Optimal Calibration Accuracy for Gravitational Wave Detectors 
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Calibration errors in the response function of a gravitational wave detector degrade its ability 
to detect and then to measure the properties of any detected signals. This paper derives the 
needed levels of calibration accuracy for each of these data-analysis tasks. The levels derived here 
are optimal in the sense that lower accuracy would result in missed detections and/or a loss of 
measurement precision, while higher accuracy would be made irrelevant by the intrinsic noise level 
of the detector. Calibration errors affect the data-analysis process in much the same way as errors 
in theoretical waveform templates. The optimal level of calibration accuracy is expressed therefore 
as a joint limit on modeling and calibration errors: increased accuracy in one reduces the accuracy 
requirement in the other. 



I. INTRODUCTION AND REVIEW 

The response function is used to convert the elec- 
tronic output of a gravitational-wave detector into the 
measured gravitational-wave signal. This response func- 
tion is determined experimentally by performing a series 
of measurements when the detector is offline, and then 
monitoring the output of the working, resonant detec- 
tor (in a time and frequency dependent way) as it re- 
acts to inputs designed to simulate its interaction with 
gravitational waves The calibration procedure pro- 
duces a response function that is known therefore only to 
the level of accuracy with which these various measure- 
ments are performed, and only to the extent the state 
of the detector changes predictably between calibration 
measurements. This paper evaluates the effects of these 
response-function errors on the subsequent gravitational- 
wave data-analysis process, and from this determines the 
optimal levels for calibration accuracy. 

Inaccuracies in the response function degrade the abil- 
ity to detect gravitational-wave signals in the noisy data 
stream; and once detected, they also reduce the ability 
to measure the physical properties of the gravitational- 
wave source that produced the signal. Errors in the 
gravitational-waveform models used in the data-analysis 
process also degrade the detection and measurement pro- 
cedures in a very similar way. An earlier discussion of 
the role of calibration error on these data-analysis func- 
tions, cf. Ref. adopted the viewpoint that the level 
of calibration error was fixed. The analysis there focused 
on determining the point at which further reduction of 
waveform-modeling errors would be made irrelevant by 
the presence of calibration error. A more proactive view- 
point is adopted here: that both the calibration error and 
the waveform-modeling error levels can (in principle) be 
set to any desired level. This paper determines the op- 
timal levels for the combined calibration and waveform- 
modeling errors needed to perform detections and also 
to perform measurements on any detected gravitational- 
wave signals. These error levels are optimal in the sense 
that lower accuracy levels would reduce the quantity and 
quality of the scientific information extracted from the 
data; while higher accuracy would be made irrelevant by 



the intrinsic noise level of the detector. 

Let us begin by discussing briefly some of what is al- 
ready known about the effects of calibration error. To 
that end, let us first establish some notation. Let v{f) 
denote the direct electronic output of the detector, and 
R{f) the response function used to convert this raw out- 
put to the inferred gravitational- wave signal h{f): 



h{f) = R{f)v{f). 



(1) 



For simplicity, the discussion here is expressed in terms 
of the frequency-domain representations of the various 
quantities. For example the frequency-domain waveform, 
h{f), is related to its time-domain analog, h{t), by the 
Fourier transform: 



Hf) 



h{t)e~^'''f'dt. 



(2) 



This transform follows the convention of the LIGO Sci- 
entific Collaboration [3] (and the signal-processing com- 
munity) by using the phase factor e~'^'^^f*, while most of 
the early gravitational- wave literature and essentially all 
other computational physics literature use e^'^'-^*. This 
choice does not affect any of the subsequent equations in 
this paper. 

Let us assume that the measured response function 
R{f) differs from the correct exact function Re{f) by 
SR{f) = R{f) — Reif)- This error in the response func- 
tion will affect measurements in two ways. The response 
of the detector to a gravitational- wave signal he will pro- 
duce an electronic output v^- So the first effect of using 
the measured response function R, is to interpret the sig- 
nal as the waveform h ~ Rve = /iee''-^"+'**" , where the 
logarithmic response function amplitude Sxb. and phase 
6^R errors are defined by 



R^ Re + SR^ ReC 



(3) 



This will produce a waveform error, 

ShR = /iee^'^«+'^*« - /le « heiSxR + iS<^R), (4) 

caused by the calibration error of the detector. The sec- 
ond effect of calibration error on measurements made 
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with the detector are errors in understanding the charac- 
teristics of the detector noise. In particular, the measured 
power spectral density of the noise Sn will differ from the 
exact Se due to the calibration error SR. The measured 
power spectral density of the noise Sn is related to 5*6 by 



(5) 



Both the detection and the measurement of a grav- 
itational wave's properties are adversely affected by 
response- function induced errors in the waveform, SHr, 
and the measured noise spectrum, S'„(/). Similar ad- 
verse effects are caused by errors in the waveform models 
used as part of the gravitational- wave data-analysis pro- 
cedure. Let 6hm{f) = hm{f) — h^if) denote the differ- 
ence between a model gravitational waveform hm, (e.g., 
one produced by a numerical-relativity simulation) and 
the exact waveform he- Both types of waveform error, 
Sh^ and 5hm, cause reductions in the signal-to- noise ra- 
tio, pm, obtained when a signal is projected onto the 
model waveform. Keeping terms through second-order 
in 5hji and 5hm, it was shown previously 0] that the re- 
sulting measured signal-to- noise ratio, pm, is related to 
the optimal signal-to-noise ratio, p, by the expression: 



Pm = p - -^{{5h,n - 5hR)j_\{Sh„ 
2p 



5h, 



(6) 



where the quantity {Shm — Shf()± is the projection of 
— 6hR orthogonal to the exact waveform, 

[ohm-OhR)^ = dh,n~OhR-he 77— m ■(') 

{rie\ne) 

The noise- weighted inner products, e.g., {ShmlShji), used 
in these expressions are defined with respect to the mea- 
sured power spectral density of the noise Snif)' 

(SKilShji) = 

Vo SM) 

The derivations of these expressions are given in some 
detail in Sec. Ill of Ref. 0]. 



II. CALIBRATION ACCURACY STANDARDS 

The expression for the difference between the measured 
and optimal signal-to-noise ratios Sp = pm — p in Eq. Q 
is remarkably simple, depending only on the difference 
between the waveform errors, 5p = Sp{Shm — Shji). At 
the most basic level, the waveform-accuracy standards 
developed in Ref. were obtained by limiting the size of 
the waveform errors to those producing acceptably small 
changes in Sp. Since Sp depends only on the difference in 
waveform errors, Shm — Sh^, from Eq. ([5]), it follows that 
the ideal-detector waveform accuracy standards can be 
extended to the realistic-detector case [Shu ^ 0) simply 



by replacing Shm with Shm — Shu in those ideal-detector 
standards. Thus the optimal accuracy requirement on 
the combined (calibration plus modeling) waveform er- 
rors that ensures no loss of scientific information during 
the measurement process is. 



(Shm - ShnlShm ~ Shji) < 1. 



(9) 



This is the generalization of the ideal-detector condition 
derived as Eq. (5) of Ref. 3. Similarly the optimal accu- 
racy requirement on the combined waveform errors that 
ensures no significant reduction in the rate of detections 
is. 



{{Shm - ShR)±_\{Shm - Shii)±) < 2p 



(10) 



The parameter emax determines the fraction of detections 
that will be missed as a result of calibration and model- 
ing errors, as discussed in some detail in Ref. i. As m 
the ideal-detector case, it is more convenient to convert 
this optimal accuracy requirement for detection into the 
slightly stronger sufficient condition. 



{Shm - ShnlShm - Sh^) < 2p 



(11) 



which does not require a knowledge of the projection 
(Shm — Shji)±. This simpler expression is the generaliza- 
tion of the ideal-detector condition derived as Eq. (15) of 
Ref. J^. 

Both of these accuracy standards, Eqs. Q and (fTTj) . 
on the combined (calibration plus modeling) waveform 
errors are conditions on the noise-weighted norm of 



Shj, 



Shfi. The waveform-modeling errors, Shm, have 



an entirely different source and are therefore completely 
uncorrelated with the calibration errors, Sh^. It is useful 
therefore to express these accuracy standards in a form 
that isolates each type of error. This can be done with a 
simple application of the Schwarz inequality, cf. Eq. (44) 
of Ref. 0: 

{Shm - ShnlShm ~ Shu) < 

(^^/{ShmlShm) + V{ShR\Shi^)y . (12) 

This inequality is fairly tight, in the sense that equality 
is actually achieved for the case Shm = —Shu- Based on 
this inequality, the following slightly stronger, sufficient, 
versions of the accuracy requirements can be constructed 
for measurement. 



a/ {Shm\Shm) + ^/ {Shii\Sh„) < 1, 



(13) 



and for detection, 



^/{Shm\Shm) + \/{ShR\ShR) < V2en,ax/0. (14) 

These conditions reduce to the model- waveform accuracy 
standards derived in Ref. Q for the ideal-detector case 
{Shu = 0). These accuracy standards place more strin- 
gent conditions, however, on the waveform-modeling er- 
ror when there is a non-negligible level of calibration er- 
ror. 
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The allowed error levels, due to calibration and wave- 
form modeling, can be apportioned between the two er- 
ror sources in any way that is consistent with Eqs. (jl3p 
and (|14p . Determining the most efRcient way to do this 
would require an analysis of the relative costs of improv- 
ing the accuracies of each error source. It seems likely 
that adopting standards with comparable requirements 
for each type of error will be close to optimal. Let us 
explore in some detail then what the resulting calibra- 
tion and modeling accuracy requirements would be in 
that case. From Eqs. and (fH|) it follows that the 
appropriate limits become 

{5h^\5h^) = {5hR\5hR) < (15) 

for measurement, and 

{Sh„,\Sh^) = {ShRlShR) < ^-^p\ (16) 

for detection. These waveform-modeling standards are 
stricter by a factor of two than those derived in Ref. Q 
for the ideal-detector case. 

It is useful to translate these accuracy requirements 
into a more familiar language, by noting that (to lowest 
order) the waveform error can be written in terms of log- 
arithmic amplitude and phase errors: 5hm ~ he{Sxm + 
iS^m). It follows that the norm of the waveform- 
modeling error can be expressed in the form, 

{Sh^\Sh^) = p2 (s^^ + W^^^ , (17) 

where the signal and noise weighted averages of the am- 
plitude and phase errors are defined by, 

Jo P^SnU) 

■5^^ = r(5$™)^4NLd/. (19) 
Jo P t>nU) 

Note that the weight term, A\he\'^ / p^Sn, which appears in 
these definitions has integral one; so these are true (sig- 
nal and noise weighted) averages of 5xm a-nd S'^m- The 
averages of the calibration amplitude and phase errors, 
are defined analogously. In terms of these averages then, 
the waveform-accuracy standards of Eqs. ((T5)) and (flBl) 
become 

v^^?vS?= \/^^?T^^?< ^, (20) 

for measurement and 

V^'+^' = \/^'+M^'< (21) 
for detection. 

For Advanced LIGO the maximum signal-to-noise ra- 
tio for a binary black- hole signal may be as large as about 



100, so the resulting accuracy requirements sufficient for 
measurement, from Eq. (|20p . for such events are 

\/^' + = \/^' + < 0.005. (22) 

Thus the averages of the frequency-domain amplitude 
and phase errors must be at about the 0.35% and the 
0.0035 radian levels respectively for measurement. If the 
Advanced LIGO search template banks are constructed 
(as in Initial LIGO) with waveform templates spaced so 
that no point in the template subspace has a mismatch 
larger than 0.03 from some element in the bank, then 
the maximum mismatch e„iax must be chosen to be 0.005 
to ensure the signal loss rate is about 10%, cf. Ref. 0]. 
In this case the resulting accuracy requirements for both 
waveform and calibration errors sufffcient for detection 
from Eq. (pij) arc 

+ = \/^' + < 0.05. (23) 

Thus the accuracy requirements for detection are an or- 
der of magnitude less stringent than those needed for 
measurement of the strongest likely sources in Advanced 
LIGO. The required averages of the frequency-domain 
amplitude and phase errors must be at about the 3.5% 
and the 0.035 radian levels respectively for detection. 

It is easy to imagine how two different sets of model 
waveforms could be designed to accomplish the two dis- 
tinct data-analysis tasks. One set could be prepared for 
use in searches of gravitational-wave signals using the 
lower accuracy standards needed for detection. And a 
second set could be prepared with the higher accuracy 
standards needed for measurement, but only in the very 
small portion of waveform parameter space where they 
are needed for measurements on previously detected sig- 
nals. At present it seems unlikely that it will be possible 
to perform detector calibrations in a way that provides 
the lower level of calibration accuracy needed for detec- 
tions at all times, and only calibrates to the higher accu- 
racy standards needed for measurements ex post facto in 
those data segments where detections have been made. 
So it seems likely that it will be necessary to calibrate 
the detectors at the level needed for measurements, e.g., 
according to the standards of Eq. (|20p . whenever data 
is collected. In this case, the accuracy standard for de- 
tections for waveform-modeling error could be relaxed 
somewhat to the level 

+ < V"!^. - ~ 0.095, (24) 

which is almost identical to the ideal-detector require- 
ment derived in Ref. 

III. USER FRIENDLY STANDARDS 

The accuracy standards derived in Sec.|TT]for the com- 
bined (calibration plus modeling) waveform errors are 
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FIG. 1: Curves illustrate C, the ratio of the standard signal- 
to-noise measure p to a non-standard measure defined in 
Eq. I|27l) . as a function of the total mass for non-spinning 
equal-mass binary black-hole waveforms. Dashed curve is 
based on the Initial LIGO noise spectrum [^; solid curve is 
based on an Advanced LIGO noise curve [^. 



not easily applied using the basic expressions given in 
Eqs. ((20|l and ([2T|l . Model waveforms are generally con- 
structed in the time domain (e.g., by performing numeri- 
cal simulations), so verifying the basic frequency-domain 
standards using estimates of the time-domain errors is 
not straightforward [3| . Neither can the basic expressions 
for the standards on the calibration errors be enforced 
in a straightforward way. While good estimates of the 
frequency-domain response-function errors are generally 
available [H, the accuracy standards of Eqs. ([20|l and ((2T|) 
require computing their averages weighted by the grav- 
itational waveform h^. What waveform should be used 
when applying these standards? This section transforms 
the basic accuracy standards of Eqs. (I13p and (fT4|) into 
forms that are more easily used by those responsible for 
calibrating the detector, and by those responsible for con- 
structing and verifying the accuracy of model waveforms 
as well. 

The norm of the waveform- modeling error {5hm\6h„i) , 
which appears in the accuracy standards of Eqs. p3|) and 
(I14p . is constructed from the frequency-domain estimates 
of those errors. It is not straightforward to determine 
useful estimates of these frequency-domain errors from 
the time-domain waveform errors that are directly acces- 
sible to the waveform-modeling community. It is useful 
therefore to transform the expression for the limits on the 
modeling error into ones based directly on time-domain 
estimates of the errors. This can be done, following the 
argument in Sec. II. C of Ref. 0, using an application of 



{Sh,n\Sh„i) < 



C2 WKitW 



(25) 



where is the norm of 5hm{t)^ defined as 

WSh^itW = / \5hM'dt, (26) 



and where C is the ratio of the standard signal-to-noise 
measure p to a non-standard measure: 



2||fee(t)|P 

min5„(/) 



(27) 



Figure [T] illustrates C for non-spinning equal- mass binary 
black-hole waveforms (cf. Fig. 4 of Ref. [2] ) . This quan- 
tity can be evaluated in a straightforward way when any 
class of model waveforms is computed. The right side of 
Eq. (|25p is always larger than the noise-weighted norm 
{5hm\5hm) that appears on the left. Sufficient conditions 
for model-waveform accuracy based on the time-domain 
norm ||(5/im(i)|| can be obtained therefore by replac- 
ing {5hm\Shm) with the right side of Eq. ((25|) wherever it 
appears in the accuracy standards of Eqs. and (|14p . 

The norm of the waveform error caused by detector 
calibration errors, {5hR\5hpi) ^ is also difficult to evalu- 
ate because it depends on the gravitational waveform /le 
in addition to the purely detector-based errors 5xr and 
(5$ij. The detector calibration-error terms in this norm 
can be isolated from the gravitational- waveform terms in 
a straightforward way using the inequality 



{5hR\5hR) 



{5x1 



Sn 



< p'^ max (Sxr 



(28) 



The right side of Eq. ([281) is very easy to evaluate, and 
is always larger than the noise-weighted norm {6hji\Shii) 
that appears on the left. Sufficient conditions for cali- 
bration accuracy based on max {Sxji + S^ji) can be ob- 
tained therefore by replacing {ShulShi^} with the right 
side of Eq. wherever it appears in the accuracy stan- 
dards of Eqs. ini) and dUl). 

In some circumstances the calibration-accuracy stan- 
dards obtained using max {Sxr + <5$|j) may be much 
stronger than necessary; for example when Sxji + 
is sharply peaked at frequencies where the detector 
noise is large. In this case it may be advantageous to 
employ a different simplification of the accuracy stan- 
dards. The detector calibration-error terms in the noise- 
weighted norm {Shii\6hji) can also be isolated from the 
gravitational-waveform terms with an application of the 
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Schwarz inequality: 



JO "5" 




(29) 



This inequality can be re-writtcn in the more compact 
form, 



2 2 



(30) 



by defining a few useful quantities. The noise- weighted 
averages of Sxb. and (5$_r are defined as 



4 4rt^ 
^XR-^df, 

4 f°° . 4n2 



where the total detector noise, n, is defined as 



df. 



(31) 
(32) 

(33) 



Note that the weight, An^ /Sn, which appears in Eqs. (pij) 
and ([32|) has integral one; so these are true (noise- 
weighted) averages of 6xr and S^r. Note the unusual 
fourth power of the averaged quantity which appears in 
these definitions. This is required because that is the 
power of the error terms which appear on the right side of 
Eq. (p9)) . Finally, the quantity C that appears in Eq. pO|) 
is the ratio of the standard signal-to-noise measure p to 
another non-standard measure: 



IP'S 11 



df 



(34) 



Figure [2] illustrates C for non-spinning equal-mass bi- 
nary black-hole waveforms using standard Initial and 
Advanced LIGO noise curves 

UM- This quantity can 
be evaluated in a straightforward way when any class of 
model waveforms is computed. It isn't completely clear 
why the Advanced LIGO version of this curve is almost 
a factor of two smaller than the Initial LIGO curve. This 
may be due to the fact that the integral of |/ie(/)|^ in 
C, Eq. ([34]) . is dominated by its low frequency contribu- 
tions where /ie(/) is largest. The Advanced LIGO noise 
curve is significantly smaller in this low-frequency range, 
so these contributions are much more significant in that 
case. 

The maximum calibration-error, max {5x\ + ^^'r) i 
and the noise-weighted averages, Sxr and S^r, which 
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FIG. 2: Curves illustrate C, the ratio of the standard signal- 
to-noise measure p to another non-standard measure defined 
in Eq. p4p . as a function of the total mass for non-spinning 
equal-mass binary black-hole waveforms. Dashed curve is 
based on the Initial LIGO noise spectrum; solid curve is based 
on an Advanced LIGO noise curve. 



appear in Eqs. and (|30p depend only on information 
that pertains to the detector itself. All of the dependence 
on the waveform he in the original norm, {5hR\5hR) , has 
been moved into the signal-to-noise ratio p and the quan- 
tity C. Thus the right sides of Eqs. ^ and ^ should 
be much easier for those performing detector calibrations 
to evaluate. The right sides of Eqs. (|28p and ([50)1 are 
always larger than the noise- weighted norm {dhulShu) 
that appears on the left. Sufficient conditions for model- 
waveform accuracy based on the maximum calibration 
error max [Sxji + S^ji) (or the noise-weighted averages 

of the calibration errors Sxr and (5$^) can be obtained 
therefore by replacing {6hii\5hji) with the right side of 
Eq. or ([50)1 wherever it appears in the accuracy stan- 
dards of Eqs. and (fH)) . The resulting accuracy stan- 
dards based on the maximum calibration error [and using 
the re- written norm of the waveform-modeling error from 
Eq. ^] become. 



1 \\shi^m 
c wKm 

for measurement, and 



max {6x\ 



1 wshinim 

C \\hS)\\ 



(35) 



(36) 



for detection. Analogous versions of the standards based 
on the noise-weighted averages of the calibration errors, 
6xr and S^r, can be obtained in a similar way. Ei- 
ther version of the standards is sufficient to guarantee the 
needed level of accuracy for gravitational- wave data anal- 
ysis. The most efficient choice for a particular detector. 
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and for a particular type of source, will be determined 

-2 2 



by whether max ((5x?{ + or {8xr + '5$_r")/C^ 

smaller. 



IS 



IV. DISCUSSION 

A new set of accuracy standards have been developed 
for the calibration and modeling errors of the waveforms 
used for gravitational- wave data analysis. The basic 
standards, Eqs. and (fH)) . are expressed most natu- 
rally in terms of the noise-weighted inner products com- 
monly used in gravitational- wave data analysis. These 
basic expressions are not very convenient for actually ap- 
plying the standards, however. So the basic standards 
have been transformed into expressions that are easier to 
apply: Eqs. ([55]) and (|36p. These new representations of 
the accuracy standards are slightly stronger, and if sat- 
isfied are sufficient to ensure the original standards are 
satisfied. 

The new accuracy standards, Eqs. ([35|l and ((36|) . pre- 
scribe a maximum for the combined calibration and mod- 
eling errors of the gravitational waveforms, not for each 
type of error separately. This means that an increased 
accuracy in one allows a somewhat weaker requirement 
on the other. Determining how to aportion the accuracy 
between the two error sources in an optimal way would 
require a careful analysis of the costs involved in reduc- 
ing the error from each source. It seems reasonable to 
expect that requiring approximately equal accuracy for 
the calibration and modeling errors will be close to op- 
timal. In this case the transformed expressions for the 
new accuracy requirements arc 



1 W^hram 

c wKm 

for measurement, and 

1 \\&K,{t)\\ 



max ((5x|f, + < 



1 



2p„ 



c wKm 



max + ^^^d < 



(37) 



(38) 



for detection; the constant pmax represents the signal- 
to-noise ratio of the strongest detected source. If the 
calibration of the instrument must be maintained at the 
level needed for accurate measurements of the strongest 
anticipated sources during the entire data collection pe- 
riod, then the accuracy requirements on the waveform- 
modeling error can be relaxed somewhat for detection: 



1 wsh^m 
c wKm 



< V2e 



1 



2p„ 



(39) 



These new accuracy standards should be applicable for 
essentially any gravitational-wave detector and any type 
of model waveform used in the data analysis process. To 
apply the standards for each particular case, the quan- 
tities Pmax, Emax, C and (pcrhaps) C must be evaluated 
for the particular family of model waveforms, using the 



noise spectrum of the particular detector. Some insight 
can be gained into what the standards will actually look 
like by examining the case of binary black-hole inspiral- 
merger-ringdown waveforms using the Advanced LIGO 
noise curve. The quantities C and C have been com- 
puted for this case using an equal-mass non-spinning bi- 
nary black-hole waveform obtained by matching together 
numerical and post-Newtonian waveforms IE, l8|] . The re- 
sults are depicted in Figs. [1] and [2] for binary systems with 
total masses in the range 4 — 4OOM0. From these graphs 
we see that C > 0.019 and C > 0.23 for these waveforms 
and the Advanced LIGO noise curve. The strongest bi- 
nary black-hole signals in Advanced LIGO are expected 
to have signal-to-noise ratios that may be as large as 
Pmax ~ 100. Assuming the template banks of model 
waveforms are constructed in the same way as those for 
Initial LIGO, the maximum mismatch compatible with a 
10% event loss rate is Cmax = 0.005. Substituting these 
values into the accuracy standards of Eqs. ([ST]) and pS]) 
results in the following calibration and waveform model- 
ing accuracy requirements for Advanced LIGO: 

53^^^ " \/^^^(^xl + S<^l) < 0.005 (40) 
for measurement, and 



53 



wheim 



^ max {5x 



5^1) 



< 0.05 (41) 



for detection. If the calibration accuracy is fixed at the 
higher level needed for measurements in the strongest 
sources for the entire period in which data is collected, 
then the standard on waveform-modeling error for detec- 
tion can be relaxed to 



53^^#< 0.095. 

Wheim 



(42) 



A somewhat troubling feature of these conditions is 
the rather large coefficient 1/C w 53 that multiplies the 
norms of 5hm{t) in these expressions. This is really 
just an artifact of the extremely long model waveform 
(containing about 1000 wave cycles) used here when eval- 
uating C 1^. The quantity C contains the norm of 
the waveform heit) in its denominator, and this norm 
becomes quite large when it is estimated using model 
waveforms hm with many wave cycles. This issue is dis- 
cussed at some length in Ref. I4| and will not be addressed 
further here, since it does not bear directly on the main 
focus of this paper: deriving the optimal levels of cali- 
bration error. 
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